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Abstract 



o. 

tN ■ Gangbo and McCann showed that optimal transportation between 

hypersurfaces generally leads to multivalued optimal maps - bivalent 

Q^ ■ when the target surface is strictly convex. In this paper we quan- 

.^^ ] tify Holder continuity of the bivalent map optimizing average distance 

• ■ squared between arbitrary measures supported on Euclidean spheres. 

■*— > 

1 Introduction 

Let X, Y be two measure spaces, /i, v two probability measures defined on 
£> ■ X and Y, respectively, and c a measurable map from X xY to [0, +00]. Let 

us denote with T the set of all the probability measures onlxf that have 
jV I marginals fi and v. The minimization problem 

n: r 

t"-- ■ minimize / c(x, y)drf(x, y) among 7 € T(/i, v) 

O ' JXxY 

o_ 

is known as Kantorovich's optimal transportation problem; c is called the 
cost function, and every probability measures in r(/i, v) is called a transfer- 
ence plan. Kantorovich's problem is meant to show how a certain mass \x 
distributed on a domain X is transported to another location (described by 
v and Y) at a minimal cost (see [32] for an exhaustive description). 

When a measurable map T supports the (unique) optimizer 7 £ T(/i, v) 
for Kantorovich's problem, T is called an optimal map. Several authors 
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treated the regularity of optimal maps when the cost function is the Eu- 
clidean squared distance; among them Caffarelli [5] [6] [7] [8] [9] [ID] . De- 
lanoe [IT], and Urbas [31] . In particular, Caffarelli showed that if the domain 
Y is convex, dfi = fdVol,dv = gdVol, where dVol denotes the Lebesgue 
measure, and the densities /, 1/g are bounded, then the optimal map is 
Holder continuous. 

Recently, in [23], Ma, Trudinger, and Wang extended regularity results 
to problems with different cost functions by introducing on the cost function 
a condition (named A3 in their paper), which implies the regularity of the 
optimal map. Their work was inspired by Wang's, and Oliker and Waltman's 
work on the problem of reflector antenna [33] [25] , which has been shown to 
be equivalent to optimal transportation of measures on the Euclidean unit 
sphere with respect to the cost function —log \x — y\ (see for example [34] 
and H8J) . 

Following Ma, Trudinger, and Wang's paper, Loeper [23], Kim, and Mc- 
Cann [19] [20j clarified the role of (A3) in the regularity theory. More pre- 
cisely, under the conditions (A0)-(A3), suitable convexity hypotheses on the 
domains, and the absolute continuity of the Lebesgue measure with respect 
to v [23], Loeper was able to prove the Holder continuity of the optimal map. 
On the other hand, Kim and McCann [20] found a covariant expression of 
(A3), named (A3s) in their paper, and extended the theory of Loeper to 
transportation problems set on a pair of smooth manifolds. Very recently, 
Liu [22] improved Loeper's Holder exponent from l/(4n — 1) to l/(2n — 1), 
and showed that the latter is sharp. 

Our paper relies on the theory of Loeper, Kim, and McCann to im- 
prove the regularity results obtained by Gangbo and McCann p3j for a 
transportation problem between boundaries of convex sets. Optimal trans- 
portation between boundaries of convex sets does not generally lead to a 
single- valued map, but rather to multivalued mappings. This is the case 
of the Kantorovich problem analysed by Gangbo and McCann, who found 
a bivalent mapping. The novelty of our paper is the quantification of the 
continuity in this setting of multi- valued mappings. 

Let O and A be two bounded, strongly convex, open sets in R n+1 , with 
Borel probability measures /i on dVt and v on <9A. We consider the Monge- 
Kantorovich problem 



d (fJ,,v):= inf / \x - y\ d>y{x,y), (1) 

where r(/i, u) denotes the set of all Borel measures on R n+1 x R n+1 having 
jjl and v as marginals. Gangbo and McCann [13] showed that two maps, 



named t + and t~, are required to support the unique optimizer 7 € T(/j,,i/). 
This means that the mass at a point x £ d£l does not always have a unique 
destination on dA, but can be split into two different destinations, t + (x) 
and t~(x). Gangbo and McCann proved that the images of a source x € dfl 
depend on x in a continuous and continuously invertible way through t + and 
t~ . In particular, t + is a homeomorphism between d£l and dA. Moreover, 
they conjectured Holder regularity for t + on d£l \ So, where 

S := {x G dn I nn(a:) ■ n A (t + (x)) = 0} 

represents a part of the 'boundary' between the region where the mass splits 
and the region where it does not. More precisely, if S2 denotes the region 
where the mass splits (bivalent region), then So contains those limit points 
of S2 at which the split images degenerate to a single image. 

The peculiarity of ([T]) is the 'hybrid' setting given by the choice of the Eu- 
clidean squared distance cost for a transportation problem set on embedded 
hyper surf aces. One of the difficulties we encountered has been to combine 
the convexity notion deriving from the Euclidean cost with the dimension 
and the pseudo-Riemannian structure of the manifolds where the measures 
are supported. Since the Hausdorff dimension of sptfi and spt^ is n rather 
than n + 1, we are not able to adapt Caffarelli's regularity theory to our 
problem; (see however |13j). Nevertheless Gangbo and McCann's conjec- 
ture about Holder continuity is reinforced by Example 2.4 of Kim-McCann 
|19j : the authors showed that the Euclidean squared distance cost, in the 
settings of ([I]), satisfies (A0), (Al), (A2) and (A3s) on 

N := {(x,y) ednx8A\ n n (x) ■ n A (y) > 0}. 

Despite this comforting result, the regularity of t + is not immediate. Loeper's 
theory needs to be adapted to our 'hybrid' setting. Moreover the target mea- 
sure with respect to t + , v\, which is the portion of mass 'transferred' by t + , 
does not inherit the hypothesis on v of having a positive lower bound on its 
density with respect to the Lebesgue surface measure. This means there are 
regions in <9£1 where the Lebesgue surface measure is not absolutely continu- 
ous with respect to v\, so one necessary hypothesis of Loeper's theory is not 
satisfied. We will treat these regions separately with a different argument. 
Our paper is organized as follows. In Section [2] we report the main result 
of Gangbo and McCann's paper |14| ; we also discuss the most important 
statement of this paper and the strategy we are going to adopt to prove 
it. We will restrict our argument to the case of spherical domains, dVl = 
dA = S n , though we believe that our regularity result can be extended to 



more general uniformly convex domains. In Section [3] we introduce and 
clarify some notation. In Section [5] we comment on some questions related 
to our problem. In Section[5]we adapt Loeper's theory to our transportation 
problem, restricting his argument to the regions of d£l where the necessary 
hypothesis on the measures holds. The regularity result on the remaining 
regions is derived in Section [6] . 

2 Preliminaries, strategy, and results 

We recall the following definitions from p3]. For a smooth convex domain 
£1, strong convexity asserts the existence of a positive lower bound for all 
principal curvatures of d£l. 

Definition 2.1. A pair of Borel measures (fJ,,v) respectively on (d£l,dA) is 
said to be suitable if 

(i) there exists e > such that fi < ^V. n [dQ and v > eHJ 1 ^^, and 

(ii) Q is strongly convex. 

If the above hypotheses are satisfied also when the roles of fi <-> v and <-> A 
are interchanged, we say the pair (fi, v) is symmetrically suitable. 

Under these hypotheses on the measures, Gangbo and McCann were able 
to prove the following optimality result. 

Theorem 2.2. Fix bounded, strictly convex domains Q,A € R n+1 with 
suitable measures fi on d£l and v on dA. Then the graphs of a pair of 
continuous maps t + : d£l — > dA and t~ : S>2 — > T2 contain the support of 
the unique minimizer 7 £ r(/u, v) for (JXJ) / 

{(x,t + (x))} x&pt ^ C spt7 C {(x,t + (x))} xedn U{(x,t-(x)))} xeS2 

(=d^n(dn xdA)). 

Here ip is Brenier's convex potential for ([TJ), £2 = dQ \ domV^, T2 C dA, 
and t + = t~ on S2 \ S2, whereas t + (x) — t~(x) 7^ is an outward normal 
for dQ, whenever x € S^- If (/-*, v) are symmetrically suitable, then t + and 
t~ [§ 2 are homeomorphisms and the minimizer 7 can be expressed by 

7 = 7i + 72, 7i = (id x i + )fl/ii, 72 = (id x t~)^ 2 , 

where \i\ := (t + )7 v\, ^2 '■= /J> — fix, and v\ := v\t%, with T£ := dA \ t~(S2) 
and the identification t~(x) := t + (x) for x € Oil. \ S^. 



Gangbo and McCann [14] introduced also a decomposition of the domain 
<9fi, which gives a geometric characterization of t + and t~. 

Proposition 2.3. Given Q, A, (/i, v) as in Theorem \2.2l we decompose d£l = 
Sq, U S\ U S2 into three disjoint sets such that 

x 6 So <^=> t + (x) = t~(x) and tiq(x) ■ n\(t + (x)) = 0. 
i£5i ^=> t + (x) = t~(x) and njj(x) • 7ia(£ + (x)) > 0. 
x £ S2 ^=> t + (x) ^ i~(s) and njj(x) • tia(£ + (x)) > 0, 

n n (x) -n A (t~(x)) < 0. 

Moreover, the sets So and So U S2 are both closed, whereas S± is open.. 

The partition of d£l will play an important role in our paper, so it is 
essential to understand the meaning of these sets. The mass lying on So U 
Si is transferred without splitting to a target set on <9A by t + , while the 
mass lying on 5*2 splits into two destinations, which are described by t + 
and t~ . For this reason we will call So the degenerate set, S\ the non- 
degenerate univalent set, and S2 the bivalent set. When the measures (//, v) 
are symmetrically suitable, an analogous decomposition of dA = T0UT1UT2 
can be introduced (see Definition 3.6 of [H]). In particular T2 is the bivalent 
set for the Kantorovich transportation problem ([1]), where (S7, z^) and (A, u) 
are exchanged, with (A, v) playing the role of the source. 

Our aim is to prove that the map t + : <9£1 — > dA is Holder continuous 
on Si and S2. One of us has been able to show that t + satisfies bi-Lipschitz 
estimates when n = 1 [28], but the argument used cannot be extended to 
higher dimensions. It heavily relies on the results of Ahmad [I]. Here we are 
developing a different strategy which works for all n > 1, when dVL, dA = S n . 
We will proceed in two steps. First we will show that t + is Holder continuous 
on the preimage S n \ (t + )~ l (T2 U To) of the set where 

vi>eH n [ d A, 

where e is the constant from Definition 12.11 This lower bound on v\ allows 
us to adapt the argument used by Kim and McCann in [20] , On (t + )~ 1 (T2), 
where the lower bound fails, the regularity of t + will be derived from the 
Holder continuity of t + on 52- In the end we will be able to obtain the 
following result. 

Theorem 2.4 (Holder continuity of multivalued maps outside the degener- 
ate set). // (fi, u) are symmetrically suitable measures on (S n ,S n ), n > 1, 
then 



t+EC^iSi) and t + eC'(5 2 



3 Notation 

The notation we are going to use is similar to that of |T3] and |19j . in 
particular we refer to Example 2.4 of [19], with <9f2 = dA = S n , c:S"xS"-> 
R, c(x,y) = \x — y\ 2 , N := {(x,y) £ S" x S n ns™(x) ■ ras n (y) > 0}, and 
N(x) := {y G S™ | (x,y) G N}. We will always use the variable x for points 
on the source domain d£l = S n , and the variable y for points on the target 
domain dA = S n . 

Let us recall the usual system of local coordinates for the points of S n 

ifi : S n D{x <E S n \xi > 0} ->• R n , <fi(x) = (xi,...,Xi-i,Xi+\,... ,x n ). 

Following this example, given x G S n and y G -/V(x) we can consider a system 
ir x of local coordinates projecting on the hyperplane perpendicular to x. In 
this way both x and y can be represented in local coordinates by means of 
the same map ir x 

x —> x, y — ^ y, 

where the capital letters stand for the image of the projection. To simplify 
the notation, given a function F : R n+1 — > R and a projection ir Xo , whenever 
x € N(x ) we will write F(X) to denote F(7r~ x (X)) = F(x). We will 
therefore write ij){X), c{X,Y) instead of Vfco 1 ^))' ^(^(I), (^{Y))- 
For example, given x G S n and y S N(x), by mean of 7r x we can write 



c(X, Y) = \X- Y\ 2 + (Vl-I^l 2 - Vl-I^l 2 ) 2 - 

In local coordinates, we use the notation Dc = (q^-, • • • , gx~ and Dc = 
(^-, . . . , -§§-) to denote the partial derivatives. The cross partial derivatives 
DDc at (x, y) G N define an unambiguous linear map from vectors at y to 
covectors at x. 

Hereafter dT-L n denotes the Hausdorff measure of dimension n, J\f p (B) 
represents the /^-neighbourhood of a set B, and [Y"o,Yi] indicate the Eu- 
clidean segment whose extreme points are Yq and Y\. 

In Section[6]we will use the expression 'angle between two vectors Z\ and 
Z2 G R n+ . The term angle refers to the arccos z \? 2 1 . 

4 Some related questions 

4.1 Relation between the convex potential ij) and the map- 
pings t + , t~ 

In a Kantorovich problem, when the optimal map exists, it is known to be the 
gradient of a convex potential. So, what is the relation between the optimal 
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mappings t + , t~, and the convex potential ip? Can we derive any regularity 
for ip from Theorem l2.4l . f ' Gangbo and McCann answered to the first question 
in Lemma 1.6 of [13] • Indeed the maps t + and t~ correspond to the outer 
and inner trace of Vip, respectively. So we can write the subdifferential of ip 
in terms of the optimal mappings: dip(x) = [t + (x) , t~ (x)] at any boundary 
point x € 80,. Moreover, in Corollary 4.4 of [14] , Gangbo and McCann 
proved that, when Q is bounded strongly convex, A is bounded strictly 
convex, and (/i, v) are suitable measures on <90, <9A, then ip is tangentially 
differ entiable along d£l. This answers the second question. From Theorem 
2.41 it follows immediately that 



i/> € CJ- 1 on 5! C S n , 
i.e. on the non-degenerate univalent set, where dtp(x) = {Vip(x)} = {t + (x)}. 



Notice that the conclusion of Theorem l2.4l does not imply ip € Cj ' 4 ™ _1 on S^- 
Indeed, ip is not differentiable in the normal direction to the sphere on S'2. 
Nevertheless, choosing the coordinates of Lemma A.l of [14] . -^- exists for 



i = 2, 3, . . . ,n + 1, and 

01b 

— — (x) = t + (x)i = t~(x)i, for i = 2, 3, . . . ,n + 1, and x € S n . 

OXi 

We conclude that the restriction of ip to S n has a derivative which is Holder 
continuous locally on Si and S^- 

4.2 The regularity of t + on So 

Presently we do not have any regularity result, except continuity from [14] . 
for t + on the degenerate set Sq. On the contrary, we will see in the state- 
ments of Theorem 15.11 and Theorem 16.11 that, on S n \ So, close to So the 
Holder constant of t + provided by our proof may become very big. More- 
over, as noticed in Example 2.4 of [19] . the nondegeneracy hypothesis (A2) 
fails on So- Therefore, we cannot apply Loeper's argument on So- On the 
other hand we believe the set So to be small. In dimension n = 1, with f2 
and A bounded strictly convex planar domains, Ahmad [I] proved that So 
consists of at most two points. 

4.3 Extending the results to more general domains 

Theorem 12.41 can be extended to the problem of transporting a measure 
on a given Euclidean sphere to a measure on any other Euclidean sphere, 



possibly with a different centre and radius. Indeed, identities (9) and (10) 
of [14] indicate how to reduce this more general problem to the case treated 
in this paper. 

Thanks to the results in Example 2.4 of [19], Theorem 15.11 can be ex- 
tended to the transportation problem where the measures (/i, v) are sup- 
ported on (dfl,dA), with f2,A C R n+1 bounded convex domains with C 2 - 
smooth boundaries. 

5 t + is Holder continuous on (t + ) -1 (Ti) C S n 

In this section we are going to adapt Kim-McCann's version of Loeper's 

argument (Appendices B,C and D of [20]) to our mapping t + , which satisfies 

(t + )$ni = v\. Thus, let us recall the regularity conditions (A0),(A1), (A2), 

and (A3s) from [19] [21] on a cost function c : S n x S" ->■ R 

(AO) (Smoothness) c £ C (N), where N has been define in Section^ 

(Al) (Twist condition) A cost c € C 1 (A r ) is called twisted if for all x € S n 

the map y — > —Dc(x,y) from N(x) C <9A to T*(S n ) is infective. 

(A2) (Non-degeneracy) A cost c £ C 2 (N) is non- degenerate if for all 

(x,y) € N the linear map DDc : T y S n — > T*S n is bijective. 

(A3s) (Strictly regular costs) A cost c E C A {N) is strictly regular on N 

if it is non degenerate and for every (x,y) € A^ 



sec. 



{x,y) {P © 0) A (0 © p) > for all null vectors p © p <E T^) N, (2) 



and equality in (0j implies p = or p = (see \19$ \2$ for an explanation 
of the notation sec; we do not need to specify it here, for it will not reappear 
in the sequel). 

Assuming (/i, u) to be suitable measures on (S n ,S n ), in order to apply 
Kim-McCann's argument we need v\ to satisfy 

there exists ei such that v\ > ci^Iqa. (3) 

From the definition of z^i in Theorem 12.21 we see that z^i satisfies ([3]) only 
outside the bivalent set T^ £ <9A = S n , i.e. outside the set where the image 
of t + is bivalent. This is the reason why we can state a regularity result 
only on a portion of the source domain, (t + ) _1 (Ti) C S n . Hereafter we will 
assume n > 1. 

Theorem 5.1. Suppose (//, v) are symmetrically suitable measures on (S n , S n ) 
(in particular, from Definition \2.1\ there exists e > such that v > eH n \_s n )- 
Then t + is locally Holder continuous on (t + )~ 1 (Ti), with Holder exponent 



at least 4n _ 1 . Our control on the local Holder constant depends on e, n, and 
tends to infinity when one approaches the boundary of N . 

Remark 5.2. Computations that show the explicit dependence of the Holder 
constant on the distance of the boundary of N can be found in \29$ . 

Lemma 5.3. The set 

N = {(x,y) G S n x S n | n S n(x) • n S n(y) > 0} 
is bi-convex in the sense of Definition 2.5 of fh 



Proof: Fix xo € S n . N(xo) appears convex from xq if and only if 
Dc(xq, N(xq)) is convex in T* (S n ). Suppose Dc(xo,yo),Dc(xo,yi) £ Dc(x,N(x)), 
where yo, y\ € N(x). We are going to show that for every 9 € (0, 1) 

9Dc(x Q , yi ) + (1 - 9)Dc(x , y ) € Dc(x ,N(x)). (4) 

Let's consider a system of local coordinates. Given xq € S n we project xq 
and y 6 N(xq) to the hyperplane perpendicular to nQ(xo) and containing 
the origin (notice that this choice of local coordinates is well defined since 
nn(x ) ■ h A (y k ) > 0, when y k <E N(x ), k = 0, 1) 

x ^0, y^Y (5) 

so that, in local coordinates, 



xo = (0,l), y = (Y,y/T=\V\ 



c(x, y) = \x- y\ 2 + (^/i-ixi 2 - Vi-iyj 2 ) 2 . 

We easily get 

Be 

-(0,Y) = -2Yi. 



dXi ^ 



If v € T x (dfl) and Vi are its coordinate with respect to the basis ^-, I can 
write 

n dc 
Dc(v)(x ,y) =v(c)(x ,y) = ^2vi——(0,Y). 



dX t 



Hence we can compute 



9Dc(v)(x , yi ) + (1 - 9)Dc(v)(x ,y ) = 6v(c)(x , y{) + (1 - 9)v(c)(x ,y ) 

= EILi [^^(o,^) + (i - e) Vi ^:(o,Y ) 

= £2=1 -2vi(9Y lti + (1 - 8)Y 0)l ). 



(6) 



Therefore, for all 6 G (0, 1) 

eDc(x Q , yi )+(l-6)Dc(x ,y ) = Dcixo^ieYi+il-O^o) € Dc(x ,N(x Q )). 

Since xq is an arbitrary point of S n , we conclude that N is vertically convex. 
By a similar argument, it is easy to show that N is also horizontally convex. 
We conclude that N is bi-convex. □ 

Proof of Theorem I5.il : Fix (x,y) = (x,t + (x)) € N, with t + (x) € T\. 
Since T\ is open, and t + is continuous, we can choose R and then r small 
enough that B r {y) C t + {Bn(x)) C T\; as asserted by Trudinger and Wang 
in [30], since N is bi-convex, taking i? and r even smaller, P = Br{x) x 
B r (y) C iV is bi-convex (alternatively, we could show directly that P is 
bi-convex, by means of the same argument used for N in Lemma 15.31 We 
replace v\ with its restriction u[ to B r (y) and we denote y! x = si(i/(). 
Taking R and then r even smaller than before gives us local coordinates 
over both domains simultaneously (for example through the chart ir x ). Let 
X = ir x (x),Y = ir x (y), and P' = ir x (B R (x)) x TV x (B r (y)). Since P is bi- 
convex and the notion of bi-convexity is coordinate invariant (as manifest 
from Definition 2.5 of [19]), P' is bi-convex with respect to the cost 



c(X, Y) = \X- Y\ 2 + (VI - |X| 2 -y/1- \Y | 2 ) 2 . (7) 

Kim and McCann showed that the cost in the original coordinates satisfies 
condition (A2) and (A3s) (see Example 2.4 of [H]), and that the quantities 
in these conditions have an intrinsic meaning independent of coordinates, 
since they are geometric quantities (i.e. pseudo-Riemannian curvatures in 
the case of (A3s) and non-degeneracy of the metric in the case (A2)). This 
implies that also the cost ((7|) satisfies (A2) and (A3s). Only the constant Cq 
of (A3s) will depend on the coordinates. Since we know that the equation 
D x c(X,Y) = Dij){X) has at most two solutions, Y + = t + (X) and Y~ = 
t~(X) and only Y + lies in P', the cost satisfies (Al) on P' . 

At this point we can apply Theorem D.l. of [20] to the cost ([7]) on P', 
with probability measures /if and vf, on tt x (Br(x)) and ir x (B r (y)) respec- 
tively, defined by 

Mi : = KOtt/^i, v{ ■■= (^)pi- 

The source fj,f is supported (and bounded above) in tv x (Br(x)) and target 
vf supported (and bounded below) in n x (B r (y)), We deduce the existence of 
a locally Holder continuous optimal map pushing jj-y forward to v[. By the 
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uniqueness of optimal transport, this map must coincide ^-a.e. with t + . 
Since both maps are continuous they agree on the (closed) support of \d x . 
Since spt/^ contains a small ball around x, this shows t + is locally Holder 
at x. □ 

6 t + is locally Holder continuous where its image 
is bivalent 

The previous section established local Holder continuity for the outer map 
t + = (s + ) _1 on the source domain s + (Ti) C S n , but not on s + (Tq U T-i) = 
So U s + (T2). Our strategy for extending this estimate to s + (T2) is de- 
scribed at the end of this paragraph. First note, however, that Gangbo 
and McCann's Sole Supplier Lemma, 2.5 of |14] , implies the outer image of 
the bivalent source is contained in the univalent target t + (S2) C T\, and 
similarly s + (T2) C Si. Since s + : S n — > S n is a homeomorphism, from 
Si U £2 = s+ (^i) U s + (T2), it follows that the bivalent source S2 C s + (Ti) 
belongs to the domain where Holder continuity of t + has already been 
shown. On this bivalent set S2, the inner map t~ is related to the outer 
map t + (x) = t~(x) + X(x)x by the geometry of the target. In Proposition 
16.21 this relation will be used to deduce (i) Holder continuity of t~ from that 
of t + . This quantifies injectivity (ii) of the inverse map s~ = (i")" 1 (through 
a bi-H61der estimate in Proposition 16. 6p . whose relation to the outer map 
s + (y) = s~(y) + uj(y)y is then used in Proposition 16. 101 to quantify injectiv- 
ity (iii) of s + = (t'^)" 1 on the bivalent target T2 = ^"(S^). This yields the 
desired local Holder continuity of t + on the source set s + (T2) mentioned at 
the outset. 

Let us recall the geometric characterization of t + and t~ from Definition 
12.31 Remembering that, on S n , ngn(x) = x, we have 

• If x G So then x ■ t + (x) = 0. 

• If x & Si then x ■ t + {x) > 0. 

• If x € S2 then x ■ t + (x) > and x ■ t~(x) < 0. 

We are going to introduce a geometric approach, based on the previous 
characterization, which allows us to prove the following theorem. Hereafter 
we will assume n > 1. 

Theorem 6.1. If (fJ,,v) are symmetrically suitable measures on (S n ,S n ), 
then t + is locally Holder continuous on (£ + )~ 1 (T2). 
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From Lemma 1.6 of [TJ] we know that t + and t are related by 

Vx € S 2 C S n t + {x) - t~(x) = \{x)x, 

where A is a continuous positive function on S 2 . Given xo,xi in S 2 we then 
have 

|t~(xi) - r(x )| < |t + (xi) - i + (x ) - A(xi)x! + A(x )x |. (8) 

We would like to exploit the regularity of t + on S 2 C (t + ) _1 (Ti), proved 
in the previous section, to prove that also t~ is Holder continuous on S 2 . 
For this purpose we also need to estimate the term \{x\)xi + \{xq)xq. This 
will be done applying the Mean Value Theorem to a suitable function and 
utilizing the geometric properties of the target. 

Proposition 6.2 (Holder continuity of t~). If t + € Cg c (S 2 ) then t~ G 
C? oc (S 2 ). Let U CS 2 andO< % := min{-x • t+(x) \xEU}. If C$ bounds 
the Holder constant for t + on U , then 



1 
ku 



C u--=(l + -r-)(C+ + 2) 



is the Holder constant for t on U . 



Proof: The function h(y) := d(y,S n ) = 1 — \y\ is differentiable on 
A = Bi(0) except at y = 0. Notice that h{t~(x)) = h(t + (x) - X(x)x) = 
h(t + (x)) = whenever x € S 2 . Consider a neighbourhood U C S 2 and the 
corresponding ku, C^ from the statement of Proposition ^. 2i Let xq, x\ G U, 
\x\ — xq\ < 2 (we need V/i to be well defined on the line segment between 
t~(xo) and t~(xi), i.e. ^ [t~(xo,t~(xi)]). Applying the Mean Value The- 
orem, we get 

= h{t + {x\) — A(xi)xi) — h(t + (xo) — A(xo)xo) 
= Vh(u) ■ (t + (xi) - t + (x ) - A(xi)xi + A(x )x ), 

for some u on the line segment between t _ (xo) and t~(x\). It follows 

(A(xi)xi - A(x )x ) • Vh(u) = (t + (xi) - t + (x )) • Vh(u). (9) 

We can rewrite ([9]) as 

(t + (xi) - t + (x )) ■ Vh(u) + A(x )(x - xi) • Vh(u) 
= (A(xi) - A(x ))xi • Vh(u); 
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then, using |V/i(n)| = 1, 

\X(x l )-X(x )\\x l -Vh(u)\ 
< |(i + (xi)-i + (xo))| + A(xo)|x -xi|. (10) 

We now state a claim, whose demonstration is postponed to the end of 
this proof. 

Lemma 6.3. Under the hypotheses of Proposition [KM, fix e € (0,1), such 
that e 2 < -4£-. Since t~ is uniformly continuous on §2, there exists 5 e , 
depending on the data through tp, such that 

\xi — xq\ < 5 e => \t~(xi) — t~(xo)\ < e. 

Then, taking xq,xi such that \x± — xq\ < 5 t , we have 

Xi ■ Vh(u) > — > for i = 1,2. 

Recalling that X(x) < 2, since d£l = S n , by means of Lemma 16.31 we 
simplify (fTOJ) to 

|A(xi)-A(x )| 

< — [\t+{ Xl ) - t + (x )\ + A(x )|x - Xl|] 
ku 

< T -[\t + {x 1 )-t + (x )\ + 2\x 1 -x \]. (11) 
ku 

Therefore, by © and (HI]), 

|£~(xi) -t~(x )\ 

< \t + ( Xl ) - t + (x )\ + A(xi)|xi - x | + |A(xi) - A(x )| 

< |t+(xi) - t + (x )| + 2|x! - x | + |A(xi) - A(x )| 

" ( 1 + ^-)l f+ ( 2; i)- t+ ( a; o)l + 2U + ^-J |xi-x |. (12) 

Combining ([12]) and t + £ C a (U), we conclude 

\t-(xi)-t-(x )\ (13) 
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i.e. t is Holder continuous on 5*2 whenever \x\ — xq\ < 5 e , with e 2 < -#■. 
We can take S e < 1, so that (fT3|) implies 

\t-( Xl )-t-(x )\ < (l + A V<7+ + 2] |xi - x | a 

= Cy\x\ — rro| Q . □ 

Proof of Lemma 1 6'. ffi Let Zj = t~(xi),i = 1,2. Notice that V/i(u) = 
— A. We have u = s^i + (1 — s)zo for some s £ (0, 1). Hence, there exists 
£ £ (0, e) such that 

X\ ■ u < —kijs + (1 — s)x\ ■ zq 
= -k v s + (1 - s)x 1 ■ (z\ +£(z - zi)) 
< -ku + (l-s)£e < -ku + e 2 . 

Using a similar argument for xq ■ u, we conclude that if e 2 < -y- then x% ■ 
Vh{u) > l^j > !%■ > 0, for i = 1, 2. □ 

Remark 6.4. Proposition \6.S\ admits a converse, i.e. if t~ £ C" oc (S2) then 
t + £ C^ (S2). This can be proved with minor changes in the preceding 
argument. 

Remark 6.5. By means of Theorem I5.il and Proposition ^.^ t~ is indeed 
locally Holder continuous on S2 with exponent An _ 1 - 

The injectivity (ii) of the inverse map s _ = (t - ) -1 on T2, is an immediate 
consequence of the local Holder continuity of t~ on S2, and it has been 
included in the following proposition. 

Proposition 6.6 (Quantifying injectivity of s~). Let V C Ti- Under the 
hypotheses of Theorem \6.1\ s~ : = (t - ) -1 satisfies 

Vyo,yi £ V sufficiently close, \s~(yi) - s~(y )\ > Cy\yi - y | 4n ~\ 

where 

Cy = (Cjj) , 

with U = s~(V) and < ky := min{— y ■ s~(y) \ y £ V}. 

Proof: Since s~ := (i")" 1 is uniformly continuous on T2, given 5 € > 
there exists 7^ > such that, if \y\ —yo\ < 7<5 e , then \s~(yi) — s~(yo)\ < 8 e . 
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Supposing 1 2/1 — j/o I < 75 s > we can apply Proposition 16 . 21 to x± = s (yi), xq = 
s~(xq) to get 

| S ~( m ) - S -(y )| > T^bi - Z/ol 4 "" 1 - □ 

We now state an elementary Lemma about vectors in R n . 

Lemma 6.7. Let u, v £ R n . Suppose the angle between u and v is less than 
\ + a, with a 6 [0, ^) . Then \u + v\ > \u\ cos a. 

Proof: Let 6 U ^ V denote the angle between u and v. Keeping |u| and \v\ 
fixed, \u + v\ can be seen as a function of 9 u ,v by mean of 

|w + f| (9 u ,v) = M + \v\ +2\u\\v\cos9 UjV , 

When 9 UtV € [0, ^ + a] , the function |-u + v\ (9 U:V ) reaches its minimum at 
8u,v = f + ct- To our purpose we can take 9 UyV = | + a. For simplicity 
we assume v parallel to e± S R n . Let us consider the projection p on the 
hyperplane perpendicular to e\ and containing the origin. Then p(u + v) = 
p(u) = \u\ cos a. Since \p(u + v)\ < \u + v\, we have the thesis. □ 

This Lemma turns out to be the key to the proof of step (hi). Under the 
hypothesis of symmetrically suitable measures, the optimal transportation 
problem we are studying is symmetric, hence every result that holds for t + 
on S n implies an analogous result for s + on S n . In particular, from Lemma 
1.6 of [H], for every y € T2 we can write 

s + (y) -s~(y) =u{y)y, (14) 

where w is a nonnegative function on T2. Hence 

\s + (yi) - s + {yo)\ = |s"(yi) - s~(y ) + uj(yi)yi - uj(y )yo\. 

If we were allowed to apply Lemma [6. 71 to the right hand side of the previous 
equality, with u = s~(yi) — s~(yo) and v = u{yi)y\ — u(yo)yo, we would 
then be able to exploit the regularity of s~ to prove step (hi). Therefore, 
we need to understand the behaviour of the angle between s~(yi) — s~(yo) 
and u(yi)yi — uj(yo)yo, when yo gets close to y\. From the monotonicity of 
dij) we have 

(s~(yi) -s~(yo)) ■ (yi -2/0) > Vyi,y e T 2 , 

which says that the angle between s~(y±) — s~(yo) and y± — yo is in [0, §] . 
If we can show that the angle between y\ — yo and oj{y\)yi — oj(yo)yo is in 
[0, a], for a certain a € [0, ^), then we can apply Lemma [6.71 to get the 
desired estimate on \s + (yi) — s + (j/o) I- 
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Lemma 6.8. Given yo,yi G Ti we denote with /3(yo,yi) the angle between 
2/1 — 2/0 and w (2/1)2/1 — ^(2/0)2/0- If the angle between yo and y\ is equal to 7 
then 

IT — 7" 



/3(2/o,2/i) € 



0. 



(15) 



Proof: The angle between yi and — yo is equal to 7r — 7, while the 
angle between y\ (or — yo) and yi — yo is ^-rp-- Since ui(yo),u(yi) > 0, 
/3(2/o,2/i)€ [O,^ 1 ). □ 

Lemma 6.9 (Dichotomy). Fix y\ G T 2 . For every integer m > 1 define 

e m (yi):=|yGT 2 |/3(y,y 1 )G \~^\ 

Unless m (yi) is empty for m sufficiently large, there exist M > and 
X > suc/t i/iat 



|s + (yi) -s + (y)| > i^|yi -y| 



Vy G 9 m (yi), witfi m> M. (16) 



Proof: We are interested in the sets m (yi) for m large, so hereafter we 
assume m > 50. Define 

< w m := inf {u(y) > | y G 6 m (yi)} 

and note zu m < w m+ i since m (yi) D @m+i(yi)- By elementary computa- 
tions, we have 



1^(2/1)2/1 - u{y)y\ cos /%, yi) 



(^(2/1)2/1 -u(y)y) ■ (2/1 -2/) 

I2/1 -2/1 



w(yi)yi 


• (2/1 


-2/) 


- u{y)y ■ (yi 


-y) 




w m j/i • 


(2/1 - 


|2/1 - 
-2/)- 


-2/1 

w(y)y • (yi - 


-y) 




ro m I2/1 


-2/I 2 


I2/1 - 

+ (m 


2/1 

m -u(y))y ■ 


(2/1 - 


2/) 



> 



> Tz m \yx-y\ 



12/1 -2/1 



Vy G m (yi) 



(17) 



where we used the definition of w m and the trivial inequality y ■ y\ < 1 to 
show that the term (ro m — uj{y))y ■ {y\ — y) is non-negative. Consider now 
the two vectors ui{yi)y\ — uj(y)y and (w(yi) — u>(y))yi, with y G T 2 . Their 
difference is parallel to yi — y, so they have the same projection on any 
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hyperplane perpendicular to y\ — y. This projection has length |w(yi)yi — 
u}(y)y\sm(3(y,yi). Therefore 

\u(yi)yi-u(y)y\sm-P(y,yi) <\u(yi)-w(y)\ VyGT 2 . (18) 

Putting together (|17|) and (|18p . we obtain an estimate for tan (~ — -) 
tan (1 - I) < tanfly,*) < M^M V y G 6 m(yi ). 

V 2 W ^m|yi-y| 

Asm-> +00, tan (f — ^) — > +00; then for every N > there exists 
m^v > 50 such that 

\u(yi) - u(y)\ > Nw m \yi - y\, Vy G 9 m (yi),m > tun- (19) 

From (|14p we have, for every y G T 2 , 

•s + (yi) - « + (y) - w(y)(yi - y) = s~(yi) - s~(y) + (w(yi) - u(y))yi. 

We define 

A := |s + (yi)-s + (y)| + |w(y)(y 1 -y)j 

> |s-(i/i)-a-(i/) + (w(i/i)-w(i/))i/i|, yGT 2 . (20) 

Using \v — u\ > |v| — |u| Vu, t> € R n+1 , we get two different estimates for A 

A > \s-( yi ) - s-(y)\ - \co( yi ) - u(y)\, (21) 

A > \ u ( yi ) - u(y)\ - |*-(i/i) - s-(y)\. (22) 

By the symmetry of the problem, using (jlip . we have 



|w(yi) - u(y)\ < ^r [\s+(yi) - s + (y)\ + 2\ Vl - y\] Vy G 9 m (yi), 

where < k' m := inf {— y • s~(y) | y G 9 m (yi)} < A^i+l- From (f2~T1) it follows 
^4 > |s"(yi) - s-(y)| - A [|s + (yi) - s+(y)| + 2|y a - y\) . (23) 



On the other hand, combining (|19p and (|22p 

-4 > Nw m \yi ~y\ - \s~(yi) -s~(y)\, Vy G 9 m (yi),m > m N . (24) 
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We can sum ([23]) and ([M]) to get 

2A > Nw m \ yi - y\ - — [|s+(yi) - s + (y)| + 2| yi - y|] 
From the definition (1201) of A, this becomes 



for every y € © m (yi), w& > m.zv- Since neither ro m nor A;^ is decreasing as a 
function of m, taking N large enough ensures N > I w^ — \- 4 ) — - — to yield 

\ K m N J ^m N 

a positive constant 

Nw mN - 4(p^ + 1) 

2 ( 1 + w L - 

\ K rn N 

such that 

\s + {yi) - s + (y)\ > K\y x - y\, Vy 6 @ m {yi),m > m N . 

To conclude we take M = mj\r- □ 

The injectivity (hi) of s + = (i + ) _1 on the bivalent target T2 = t~(S2) 
follows from Lemma 16.71 and Lemma 16.91 



Proposition 6.10 (Quantifying injectivity of s + on the bivalent target). Let 
y\ € V C Ti- Under the hypotheses of Theorem \6.1[ there exists Cy > 0, 
depending on Cy, ky (from Proposition \6.6\) . and 9(y\) (from Lemma \6.9\) , 
such that, when yo is sufficiently close to y\, 



|s + (yi) - s + (yo)| > C$\yi - y 



l/i-i 



Proof: When y^ E Q(yi,9), with 9 > 9(y±) we apply Lemma l6.9l and we 
are done. Otherwise the angle between s~ (y±) — s~ (yo) and to(yi)yi—io(yo)yo 
is smaller than ^ + 6{y\ ) . Applying Lemma 16.71 we obtain 

\s + {yi)-s + {y )\ 

= I s " (2/1) - s~(yo) +u(yi)yi -u(yo)yo\ 
> \s~ (y 1) - s~(y )\ cos 9 (yi). 

Taking yo,yi sufficiently close (\y% — yo\ < js e , from the proof of Proposition 
16.61) . Proposition 16.61 implies 

\s + (yi)-s + (yo)\ >cos9{y 1 )Cy\y 1 -y \ in -\ □ 
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Proof of Theorem \6.1[ Define yi := t + (x{) € V C T^. If yo £ 0(yi>^)) 
with # > 9(yx), we have 

|yi -y \ < K\xi -x \. 

Otherwise, by the uniform continuity of t + , taking xq sufficiently close to 
x\, we have \t + (x\) — t + (xo)| < 7<5 £ and we can apply Proposition 16.101 to 
yi = t + (xi), i = 1,0 to conclude 

1 , , i 

m -2/o I < TqrFi -xol 4 ™- 1 . □ 
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